A control volume finite element method for three-dimensional three-phase flows

A novel control volume finite element method with adaptive anisotropic unstructured meshes is presented for three-dimensional three-phase flows with interfacial tension. The numerical framework consists of a mixed control volume and finite element formulation with a new P₁DG-P₂ elements (linear discontinuous velocity between elements and quadratic continuous pressure between elements). A “volume of fluid” type method is used for the interface capturing, which is based on compressive control volume advection and second-order finite element methods. A force-balanced continuum surface force model is employed for the interfacial tension on unstructured meshes. The interfacial tension coefficient decomposition method is also used to deal with interfacial tension pairings between different phases. Numerical examples of benchmark tests and the dynamics of three-dimensional three-phase rising bubble, and droplet impact are presented. The results are compared with the analytical solutions and previously published experimental data, demonstrating the capability of the present method.     

Stochastic sensitivity analysis of coupled acoustic-structural systems under non-stationary random excitations

This paper presents a new sensitivity analysis method for coupled acoustic–structural systems subjected to non-stationary random excitations. The integral of the response power spectrum density (PSD) of the coupled system is taken as the objective function. The thickness of each structural element is used as a design variable. A time-domain algorithm integrating the pseudo excitation method (PEM), direct differentiation method (DDM) and high precision direct (HPD) integration method is proposed for the sensitivity analysis of the objective function with respect to design variables. Firstly, the PEM is adopted to transform the sensitivity analysis under non-stationary random excitations into the sensitivity analysis under pseudo transient excitations. Then, the sensitivity analysis equation of the coupled system under pseudo transient excitations is derived based on the DDM. Moreover, the HPD integration method is used to efficiently solve the sensitivity analysis equation under pseudo transient excitations in a reduced-order modal space. Numerical examples are presented to demonstrate the validity of the proposed method.     

A proper orthogonal decomposition variational multiscale meshless interpolating element-free Galerkin method for incompressible magnetohydrodynamics flow

In the recent decade, the meshless methods have been handled for solving most of PDEs due to easiness of the meshless methods. One of the popular meshless methods is the element-free Galerkin (EFG) method that was first proposed for solving some problems in the solid mechanics. The test and trial functions of the EFG are based on the special basis. Recently, some modifications have been developed to improve the EFG method. One of these improvements is the variational multiscale EFG procedure. In the current article, the shape functions of interpolation moving least squares approximation have been applied to the variational multiscale EFG technique for solving the incompressible magnetohydrodynamics flow. In order to reduce the elapsed CPU time of simulation, we employ a reduced-order model based on the proper orthogonal decomposition technique. The current combination can be referred to as the reduced-order variational multiscale EFG technique. To illustrate the reduction in CPU time used as well as the efficiency of the proposed method, we applied it for the two-dimensional cases.     

A new database and benchmark of the bond energies of noble-gas-containing molecules

We have developed a new database of structures and bond energies of 59 noble-gas-containing molecules. The structures were calculated by CCSD(T)/aug-cc-pVTZ methods and the bond energies were obtained using the CCSD(T)/complete basis set method. Many wavefunction-based and density functional theory methods have been benchmarked against the 59 accurate bond energies. Our results show that the MPW1B95, B2GP-PLYP, and DSD-BLYP functionals with the aug-cc-pVTZ basis set excel in predicting the bond energies of noble-gas molecules with mean unsigned errors (MUEs) of 2.0 to 2.1 kcal/mol. When combinations of Dunning's basis sets are used, the MPW1B95, B2GP-PLYP, DSD-BLYP, and BMK functionals give significantly lower MUEs of 1.6 to 1.9 kcal/mol. Doubly hybrid methods using B2GP-PLYP and DSD-BLYP functionals and MP2 calculation also provide satisfactory accuracy with MUEs of 1.4 to 1.5 kcal/mol. If the Ng bond energies and the total atomization energies of a group of 109 main-group molecules are considered at the same time, the MPW1B95/aug-cc-pVTZ single-level method (MUE = 2.7 kcal/mol) and the B2GP-PLYP and DSD-PLYP functionals with combinations of basis sets or using the doubly hybrid method (MUEs = 1.9-2.2 kcal/mol) give the overall best result.     

Multi-variable special polynomials using an operator ordering method

Using an operator ordering method for some commutative superposition operators, we introduce two new multi-variable special polynomials and their generating functions, and present some new operator identities and integral formulas involving the two special polynomials. Instead of calculating complicated partial differential, we use the special polynomials and their generating functions to concisely address the normalization, photocount distributions and Wigner distributions of several quantum states that can be realized physically, the results of which provide real convenience for further investigating the properties and applications of these states.     

A new non-perturbative approach in quantum mechanics for time-independent Schr?dinger equations

A new non-perturbative approach is proposed to solve time-independent Schr?dinger equations in quantum mechanics.It is based on the homotopy analysis method(HAM)that was developed by the author in 1992 for highly nonlinear equations and has been widely applied in many fields.Unlike perturbative methods,this HAM-based approach has nothing to do with small/large physical parameters.Besides,convergent series solution can be obtained even if the disturbance is far from the known status.A nonlinear harmonic oscillator is used as an example to illustrate the validity of this approach for disturbances that might be one thousand times larger than the possible superior limit of the perturbative approach.This HAM-based approach could provide us rigorous theoretical results in quantum mechanics,which can be directly compared with experimental data.Obviously,this is of great benefit not only for improving the accuracy of experimental measurements but also for validating physical theories.     

Flow physics of normal and abnormal bioprosthetic aortic valves

Flow physics of transvalvular flows in the aorta with bioprosthetic valves are investigated using computational modelling. For the efficient simulations of flow-structure-interaction in transvalvular flows, a simplified, reduced degree of freedom valve model is employed with a sharp interface immersed boundary based incompressible flow solver. Simulations are performed for normal as well as abnormal valves with reduced leaflet motion that models the effect of early leaflet thrombosis. The structure of the aortic jet and the hemodynamic stresses on the aortic wall are analysed to understand the hemodynamic impacts and possible long-term clinical implications of sub-clinical, reduced leaflet motion. The simulation results have shown that the reduced leaflet motion tilts the direction of aortic jet and generates stronger flow separation and re-attachment on the aortic wall downstream from the reduced motion leaflets. The modified flow pattern increases the wall pressure fluctuation and average wall shear stress on the downstream aortic wall, and results in the asymmetric oscillatory shear index distributions, which may have long-term clinical implications such as aortic wall damage and remodelling.     

Investigation on the effect of density ratio on the convergence behavior of partitioned method for fluid–structure interaction simulation

In order to investigate the effect of density ratio of fluid and solid on the convergence behavior of partitioned FSI algorithm, three strong-coupling partitioned algorithms (fixed-point method with a constant under-relaxation parameter, Aitken’s method and Quasi-Newton inverse least squares (QN-ILS) method) have been considered in the context of finite element method. We have employed the incompressible Navier–Stokes equations for a Newtonian fluid domain and the total Lagrangian formulation for a non-linear motion of solid domain. Linear-elastic (hyper-elastic) model has been employed for solid material with small (large) deformation. A pulsatile inlet-flow interacting with a 2D circular channel of linear-elastic material and a pressure wave propagation in a 3D flexible vessel have been simulated. Both linear-elastic and hyper-elastic (Mooney–Rivlin) models have been adopted for the 3D flexible vessel. From the present numerical experiments, we have found that QN-ILS outperforms the others leading to a robust convergence regardless of the density ratio for both linear-elastic and hyper-elastic models. On the other hand, the performances of the fixed-point method with a constant under-relaxation parameter and the Aitken’s method depend strongly on the density ratio, relaxation parameter selected for coupling iteration, and degree of deformation. Although the QN-ILS of this work is still slower than a monolithic method for serial computation, it has an advantage of easier parallelization due to the modularity of the partitioned FSI algorithm.     

均匀磁场下碳离子笔形束的剂量变化分析及位置修正方法

分析外加均匀磁场对于碳离子笔形束剂量分布的影响,并考虑修正这种影响,为磁共振成像引导碳离子放射治疗的临床应用提供指导。本文利用蒙特卡罗方法模拟计算了不同能量碳离子笔形束在不同强度磁场下的剂量分布情况,发现垂直于碳离子束入射方向的均匀磁场对于碳离子笔形束射程缩短的影响很小,磁场对碳离子束的主要影响是引起束流横向偏转,特别是碳离子束布拉格峰位置的横向侧移。横向侧移程度与碳离子束的能量和磁场强度相关,根据模拟结果,得到了一个计算碳离子束布拉格峰在磁场中相对横向偏转的方程,并提出一种校正外加磁场引起的碳离子束布拉格峰横移的角度修正方法。这些结果可用于指导磁共振图像引导碳离子放射治疗计划系统的研发。